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(^ . Abstract 

f~^ ' A set of fully relativistic one-pion-excliange electromagnetic operators is developed for use in one-particle 

emission reactions induced by electrons. To preserve the gauge invariance of the theory additional pionic 

^p ' correlation operators are required beyond the usual meson-exchange current operators. Of these, in the 

r^ . present work emphasis is placed on the self-energy current which is infinite and needs to be renormalized. 

The renormalized current is expanded to first order to obtain a genuine one-pion-exchange contribution to 

7^ . the inclusive responses. 
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1 Introduction 

In a recent paper ^ we investigated the role played by pions in inclusive electron scattering from 
nuclei within the context of one-particle one-hole (Ip-lh) excitations. There we extended our 
previous work ^, ^ in which we had pursued a systematic study of relativistic effects in the nuclear 
electromagnetic responses for various kinematical regions incorporating both meson-exchange and 
isobar currents. The resulting two-body current is a consistent first-order operator, since it contains 
all Feynman diagrams built from nucleons and pions with one exchanged pion and with a photon 
attached to all possible lines; importantly, the result is gauge invariant, as explicitly proven in ^. 

In addition to the usual contact and pion-in-flight meson-exchange current (MEC) operators, 
this fully relativistic operator includes the so-called correlation operator. Usually this is not in- 
cluded in model calculations because it gives rise to contributions which are, at least in part, 
already accounted for in the starting nuclear wave functions. However, our model is based on 
an uncorrelated relativistic Fermi gas where the ground state is a Slater determinant constructed 
with (Dirac) plane waves. Within a perturbative approach we are free to consider the one-pion 
correlation contribution to the responses as arising either explicitly in the wave functions or from 
a current operator acting on unperturbed states — our approach is the latter. When the whole 
perturbative expansion is summed up one must recover the results obtained starting with an exact 
"correlated" wave function. 

In this paper we enter in more depth into analysis of the ideas in |j|], focusing on the impact 
of pions on the nuclear electromagnetic response in the Ip-lh channel, since this gives the main 
contribution for quasielastic conditions 0]-[P. As in the MEC case, the two-body correlation 
current also contributes in this channel and thus is required to obtain a consistent, gauge invariant, 
one-pion-exchange current. When acting on the ground state with the two-body operator, the Ip- 
lh matrix element associated with the two-body currents is obtained by integrating one of the two 
particle states over the Fermi sea. In so-doing two kinds of contributions are obtained. The first 
one, sometimes referred to as a vertex correction [^, arises from the pionic correlations between the 
particle and hole; the second relates to the Fock self-energy (s.e.) 0, ^] and dresses the particle and 
hole lines. This second contribution diverges, since it corresponds to a s.e. insertion on an external 
line, which, according to field theory [^, |l^, should not be included in a perturbative expansion. 
Instead one should apply a renormalization procedure to dress the external lines by summing the 
entire perturbative series of self-energy insertions. In the nuclear case this procedure leads to the 
relativistic Hartree-Fock (HF) approach. 

In other relativistic calculations |Q, ^ this contribution has been treated by introducing from 
the start a Hartree-Fock propagator in the medium, which accounts for the s.e. diagrams. In- 
medium form factors for the Ip-lh current are also defined, however neglecting any momentum 
dependence in the self-energy and effective mass. In [M] the self-consistent Hartree mean field was 
inserted into the single-particle propagator, automatically including the Pauli blocking oi NN pairs, 
whose contribution was thus included in the RPA responses. A similar treatment including a semi- 
phenomenological nucleon self-energy in the medium at the non-relativistic level was implemented 
inref. 0. 

In Q the difficulty posed by the s.e. insertion in first order was avoided by computing the 
self-energy response as the imaginary part of the corresponding polarization propagator with one- 
pion-exchange s.e. insertions on the particle and hole lines. In this way a finite result is already 
obtained in first order (one pionic line) without having to resort to the HF approach. The question 
then arises whether it is possible to obtain the same result for the self-energy response function 
starting with finite well-defined matrix elements of the current operator. 



In this paper our aim is to answer this question by constructing a renormahzed self-energy 
current corresponding to one-pion-exchange. This current acts over free Dirac spinors and leads 
to the same response functions as those obtained by taking the imaginary part of the polarization 
propagator computed to first order. It should be clear that here the concept of renormalization has a 
many-body significance, namely it amounts to a relativistic HF approximation and ignores (see ^) 
the additional vacuum renormalization due to the change of the negative-energy sea induced by 
the nuclear medium Q. 

The new current is obtained by renormalizing spinors and energies and expanding the resulting 
in-medium one-body current to first order in the square of the pion-nucleon coupling constant. 
The renormahzed quantities should be obtained by solving a set of self-consistent relativistic HF 
equations numerically. However, to first order in the expansion, it is possible to write down analytic 
expressions for the solutions and the corresponding corrections to the bare single-nucleon current 
operator can be expressed in terms of a simple electromagnetic operator. This operator accounts 
for two main effects introduced by the interaction of the nucleon with the medium: the first is the 
enhancement of the lower-components of Dirac spinors; the second is a global renormalization of 
the spinors in the nuclear medium. These effects are genuine relativistic corrections that are absent 
within an essentially non-relativistic approach Q. Actually a third type of renormalization effect 
arises from the modification of the energy-momentum relation for a nucleon, treated in first order 
of the square of the pion-nucleon coupling constant (in other approaches this effect is embedded in 
a constant effective mass [0]). 

Using the renormahzed s.e. current operator together with the MEC and the vertex exchange 
operator we prove that the full current is gauge invariant if account is taken of the change in 
energy arising from the HF renormalization to first order. In this way the results for the inclusive 
response functions agree completely with the ones obtained in ||l| with the polarization propagator 
technique. 

The present paper is organized as follows: In Section 2 we briefly revisit the full set of Ip-lh 
current operators with one pion-exchange line which contribute to the electro-excitation process, 
paying special attention to the s.e. contribution. We show the necessity of re-defining the otherwise 
infinite self-energy diagrams. In Section 3 we perform the renormalization program, thus obtaining 
the relativistic HF equations. In Section 4 we expand the renormahzed spinors and energies to first 
order in the pion coupling constant squared obtaining a new self-energy current. We also prove 
the gauge invariance of the whole set of current operators to first order. In Section 5 we compute 
the interference response of the s.e. current with the one-body current and prove its identity with 
the one derived in the polarization propagator framework. We present numerical results for the 
response functions focusing on the contributions provided by the different pieces of the self-energy 
current. Finally, in Section 6 we summarize our results and draw our conclusions. 

2 MEC and Correlation currents 

The linked, two-body Feynman diagrams that contribute to electron scattering with one pion- 
exchange are shown in Fig. 1. The first three correspond to the usual MEC: diagrams (a), (b) 
refer to the contact or seagull current, diagram (c) to the pion-in- flight current. The four diagrams 
(d)-(g) represent the so-called correlation current and are usually not treated as genuine MEC, but 
as correlation corrections to the nuclear wave function. However, again we note that our approach 
puts all correlation effects in the current operator and uses an uncorrelated wave function for the 
initial and final nuclear states. 



In this work we use Bjorken & Drell conventions [^ and a pseudo- vector -kNN coupling, namely 

/ _ 

Ti-jrNN = i'lbl^{d^(t)a)Tail , (1) 

where ijj is the nucleon field, (pa is the isovector pion field, / represents the ttNN coupling constant 
and m,r is the pion mass. The electromagnetic currents corresponding to diagrams (a)-(g) are 
obtained by computing the S-matrix element 

^fi ~ '5'/i(-Pl' -^2' -Pli -^2) — Sfi{P^,P2,P2, Pi) (2) 

for the absorption of a virtual photon by a system of two nucleons, namely for the process 

^ + Ni + N2^N[ + N^, (3) 

with Pi, P2 {Pi, P2) being the initial (final) four-momenta of the two nucleons involved. The 
electromagnetic current is then defined according to 

Sfi{P{,P^,Pi,P2) = -ieA^{Q){27T)^6{Pi + P^ - P^ - P^ - Q)i'^(P{, P^, Pi, P2) , (4) 

where Afj^{Q) is related to the matrix element of the electromagnetic field between the incident 
photon with momentum Q and the vacuum state, namely 

{0\A,iXMQ)) = A^iQ)e-'Q-^ . (5) 

In the above j^(Pl, P2, Pi, P2) is related to the matrix element of the current as follows 

{PiPi\rmPlP2) = {27Tf6^{p[ + p'2 - q - pi - p2)i^(Pi',P^,Pi,P2) . (6) 

The general relativistic expressions for the seagull (diagrams a,b), pion-in-flight (c) and cor- 
relation (d-g) currents (j^, j^, j^^j,) are reported in [Q]. There we have also proven that when 
the seagull and pion-in flight currents are multiplied by the same electromagnetic form factor Pi , 
gauge invariance is fullfilled, i.e., Qfj.{j^ +ip +irar) = 0- 

In this work we deal with the case of one-particle emission induced by the two-body currents 
mentioned above. The matrix element of a two-body operator between the Fermi gas ground state 
and a Ip-lh excited state has the form 

{ph''\fiQ)\F) ^ (^27,fs'{ci + h-p)y^=fip,h) 

= E [(pk\J^iQ)\hk) - {pk\j''iQ)\kh)] , (7) 

k<F 

where the summation runs over all occupied levels in the ground state, and thus includes a sum 
over spin and isospin and an integral over the momentum k. Note that the current matrix element 
j^(p,h.) does not contain the factor . "^ , which should accordingly be included into the phase 

space when computing the response functions. 

It is well-known (see, e.g., 0, |2|) that the direct term in eq. (|^) vanishes for the MEC and pionic 
correlation currents upon summation over the occupied states. Only the exchange term contributes 
to the p-h matrix elements. The associated many-body Feynman diagrams are displayed in Fig. 2 
and the corresponding fully relativistic expressions for the seagull, pion-in-flight and correlation 
currents are given in |1[. 
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Figure 1: Feynman diagrams contributing to the two-body current with one pion-exchange. The 
wide hne in the correlation diagrams (d)-(g) means a fully relativistic Dirac propagator for the 
nucleon. 



In what follows we restrict our attention to the correlation current given by diagrams (d)-(g) 
of Fig. 2. Following the arguments presented in |]l|, we split the total correlation current into a 
i^ertex correlation current {v.c, also referred to as "exchange" in [^, diagrams (d) and (e)) and a 
se//- energy correlation current (s.e., diagrams (f) and (g)), according to 



Jcor Js.e. ' Jv.c. 



(8) 



The vertex and self-energy currents then read 
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Figure 2: Many-body Feyninan diagrams contributing to the one-body current with one pion- 
exchange. The thick Une in the correlation diagrams (d)-(g) corresponds to a fully relativistic 
Dirac propagator for the nucleon. Diagrams (d)-(e) represent the vertex current, while diagrams 
(f) and (g) represent the self-energy current of the hole and of the particle, respectively. 



(10) 

respectively Q. 

It is important to point out that in eq. ( |lO|) the particle {p) and hole (h) are described in the 
Fermi gas by unperturbed plane waves, i.e., they are on-shell, and hence the propagators Sf{P) 
and Sf{H) diverge. On the one hand, this problem cannot be cured at the level of the current 
matrix elements. On the other, we know that the full set of Feynman diagrams (Fig. 2) is needed 
in order to get a consistent, gauge-invariant, one-pion-exchange current and hence their presence is 
required. Indeed, the s.e. diagrams yield a far from negligible contribution to the nuclear response 
functions as shown in Q , where they are computed through the imaginary part of the first-order 
polarization propagator. Hence the necessity of having to renormalize the expression in eq. ([lO| ) 
arises. 

This divergence of the diagrams (f)-(g) is reminiscent of the well-known infinity occurring in 
standard perturbative quantum field theory, when self-energy insertions in the external legs are 
included in Feynman diagrams. As is well-known, there one should renormalize the theory by 
dressing the external legs, propagators and vertices. In the nuclear matter case we assume that the 



^Note that here, although the global factor 



V^yEpEh 
associated with the internal line, has to be retained inside the sum 



has been extracted from the current, the factor 
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Figure 3: Diagrammatical representation of the self-energy current for a hole (a) and a particle (b). 
The self-energy is defined to first order as the Fock insertion shown in (c) with one pion-exchange. 



particle-physics effects are already accounted for by the physical masses and electromagnetic form 
factors. However, an additional nuclear physics renormalization that arises from the interaction of 
a nucleon with the nuclear medium should be included at the one-pion-exchange level to account 
for the self-energy diagram. 

The self-energy current in eq. ([lO|) can be written in the following form: 



Js.e. 



.(p,h) = n(p)S(P)5^(P)r^(g)n(h) +n(p)r^(Q)5^(i7)S(ii-)n(h) 
where S(P) is the nucleon self-energy matrix that in first order reads 



S(P) 



f 



Vm'i 
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-Tal^iP- ^) • 



(11) 



(12) 



This is diagrammatically displayed in Fig. 3. The s.e. matrix, shown in Fig. 3(c), corresponds to 
the Fock term of the mean-field potential (the Hartree or direct term is zero for pion exchange). 

Performing the sum over the spin (s) and isospin (t) indices, implicitly included in eq. (|l^), and 
using the commutation properties of the gamma matrices to eliminate 75, the self-energy can be 
cast in the form 

A J_(f-A-)(A--„.)(f-A-) 



S(P) = -3Ar 
mi 



(2^)= 



2^k 



[P-Kf 



mt 



(13) 



where the sum over k has been converted into an integral. Note that the self-energies S(P) and 
S(i?) appearing in eq. ( pT|) are evaluated for free particles and holes, i.e., for P^ and H^^ on-shell. 



Hence the self-energy contributions are divergent, since so are the free propagators S{P) and S{H) 
in eq. (|ll|). Therefore they should not be computed using eq. (pi]), but rather one should first 
renormalize the wave function and the propagator of the particles in the medium. This is achieved 
through the summation of the full series of diagrams with repeated self-energy insertions displayed 
in Fig. 3. 

Now the energy of a particle in nuclear matter is modified by the interaction with the medium 
and, as well, through its energy-momentum relation. Thus the associated momentum is no longer 
on-shell and therefore in the next section we shall evaluate the self-energy for off-shell particles. 
In the first iteration, corresponding to one pion-exchange, the particle P^ is off-shell, but the 
intermediate interacting hole K^ is not ^. In this case, with the help of Dirac spinology, one writes 

(f- ^)(^ - m)[f>- f{) = 2{P-K - m2)(f + m) - (P^ _ m'^){f{ + m) , (14) 

which allows one to recast the self-energy in eq. ([l^ ) for the off-shell momentum P in the form 



T(P\ ^/' f d^k 1 2{P-K- m'){P + m) - (P^ - m'){fi + m) 

Note that the second term inside the integral vanishes for P on-shell. 

In general the self-energy of a nucleon in nuclear matter can be written in the form |^]: 

S(P) = mA{P) + B{P)-fop° - C{P)-f ■ p . (16) 

In contrast to the quantum field-theory decomposition S(P) = mA + B /P, owing to the non- 
invariance under a boost of the step function 6{kF — k) appearing in the self-energy, in nuclear 
matter B[P) ^ C{P). This in turn reflects the existence of a privileged system, namely the lab 
system where the Fermi gas has pfg = 0, this being the natural one in which to compute the 
self-energy. Under a boost, the Fermi gas ground state is no longer characterized hj k < kp and 
also the self-energy takes a different form. 

In the case of the Fock self-energy in eq. ( [l5|) the functions A,B,C can be expressed in terms 
of the integrals 

r d^k ^ 1 

and in terms of the functions Ko{P) and Kz{P), defined as | 

Lo(P) = K^{P)po (19) 

Up) = K3{P)p. (20) 



Indeed one gets 



A{P) = -3i^[2(P^L'^(P)-m2/(P))-(p2-m2)/(P)l (21) 



m^ 



^Note that in deriving eq. (|l3| ) we have assumed free spinors u(k); hence eq. (|l^) is only valid for K'^ on-shell. 
The off-shell case requires one to redefine the spinors ii(k) according to an interacting Dirac equation, as is shown 
later. 

■^ L is parallel to p since, choosing p along the z axis, the azimuthal integration in eq. (hsl) yields L^ — Ly — 0. 
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Figure 4: Diagrammatical series for the nucleon propagator in the medium. 

B{P) = -2,^\2{P^L^'{P)-mH{P))-{P^-n?)KQ{P)] (22) 

C{P) = -3^\2{P^L''{P)-m^I{P))-{P^-m^)Ks{P)\ . (23) 

Note that A = B = C for P on-shell. In this case one has 

^PU-sheii = AiP){m+ n . (24) 

3 Hartree-Fock renormalization in nuclear matter 

3.1 Nucleon Propagator 

The nucleon propagator in the nuclear medium results from summing up a series with any number 
of self-energy insertions (see Fig. 4) , namely 



^ +^;r^S(P)— ^ + — ^S(P)— ^S(P)^^^ + --- 



P — m ' f — m ^ ' P — m ' f — m ^ ' f — ni ^ ' f — m ' f — m — S(P) 

Using the spin decomposition of the self-energy in eq. ([T^), we can write 

p-m- S(P) = [1 - B{P)\ 70P0 - [1 - C{P)\ 7 • p - [1 + A{P)\ m . 
Now the new four-momentum /^ = f^{P), which is related to P^ as follows 

1 - B{P) 



and the functions 



fo{P) 
f(P) 

m{P) 

z{P) 



1 - C{P) 
P , 



Po 



l + A{P) 
1 - CiP) 

1 
1 - C{P) 



m 



■ (25) 

(26) 

(27) 
(28) 

(29) 
(30) 



allow one to recast the nucleon propagator in the form 

1 z(P) 



z{P) 



p-m-^P) joMP) - 1 ■ P - m{P) i{P)-m{P) 



(31) 



For a nucleon with a fixed three-momentum p, the pole of the propagator in eq. (^) in the variable 
Pq defines the new energy of the nucleon in the medium. To find the latter we introduce 



E{P) = E{p, m{P)) = A/p2 + m(P)2 . 



Then the propagator reads 



1 



z{P)^ fjP) + mjP) 

p-m- S(P) /o(P) - E{P) UP) + E{P) 

and its pole po is found by solving the implicit equation 

UP) = E{P) , 

which, exploiting eq. (p7|), can be recast as follows 



Po 



1 - CjP) 
l-B{P) 



p2 + m(P)2 



l-g(po,p) 
l-Bipo,p) 



p2 + [m(po,p)]' 



(32) 



(33) 



(34) 



(35) 



The solution of eq. (|35| ) for fixed p defines the new dispersion relation pQ = e(p) for interact- 
ing nuclear matter. Once the above equation has been solved, the field strength renormalization 
constant 

z{P) 



Z2{p) = Res 



UP)-E{P) 



(36) 



PO=<:(p) 



defined as the residue of the first factor on the right-hand side of eq. (|33|) at pq = e(p), can be 
computed. Indeed using eq. (|30[), Z2{p) is obtained by expanding the denominator around the pole 

e(p), i.e., 

[1 - C{P)] \UP) - E{P)] = Z2{p)-' [po - e{p)] + • • • ; 



(37) 



hence 



Z2{P)-' 



d 
dpo 

d 
dpo 



PO=<:(p) 



[[l-C{P)][UP)-E{P)\} 
{[l-B{P)]po-[l-C{P)]E{P)] 



PO=e(p) 

With the help of eq. (p^) the derivative can be easily evaluated, the result being 



(38) 



^2(P)"' 



dB m dA p^ dC 

dpo E opo E opo 



(39) 



PO=e{p) 



3.2 Nucleon Spinors 



The self-energy modifies not only the propagator and the energy-momentum relation of a nucleon, 
but, as well, the free Dirac spinors. In fact the spinors are now solutions of the Dirac equation in 
the nuclear medium |jl2|, i.e., 

[F-m- I](P)]0(p) = , (40) 



which, again using the decomposition in eq. (jT^), can be recast as follows 

bofo{P)-l-p-rn{P)]cl){p)=0, 



(41) 



the functions fo{P) and m{P) being defined in eqs. (^) and (p9|), respectively. Equation (^Tj) has 
the same structure as the free Dirac equation; hence for the positive-energy eigenvalue one has 



/2(P)=p2 + m2(p) 



(42) 



which implicitly yields the energy pQ = e(p) of the nucleon in the nuclear medium. This result was 
already obtained as the pole of the nucleon propagator. Then the corresponding positive-energy 
spinors {s being the spin index) read 



^s{P) 



Z2{P) 



' E{p)+m{p) 
2m(p) 



1/2 



Xs 
crp 



^(p)+m(p) 



Xs 



Z2{p)us{p,m{p)) (43) 



and the functions of the three-momentum p 

m(p) =m(e(p),p) 



E{p)^E{eip),p) = ^p^ + mip)^ 
represent the nucleon's effective mass and effective energy corresponding to po 



strength renormalization constant, y/Z2{p), of the new spinors, defined in eq. 
renormalization theory, since the propagator in eq. ( pll) for po close to the pole e(p) reads 



(44) 

(45) 

e(p). The field 
is required by 



m{p) 



r 



m 



S(P) po - e(p) E{p) 



l]0s(p)0s(p) 



Z2{p) /(p) + m(p) 
Po - e(p) 2£'(p) 



(46) 



Once the new spinors have been computed, the self-energy can be evaluated by inserting the latter 
into eq. ([l^). Then the Dirac equation should be solved again with the new self-energy and so on. 
This self-consistent procedure leads to the relativistic Hartree-Fock model which has to be dealt 
with numerically. 

In this paper we do not attempt to solve the HF equations, since we are interested only in 
the first-order correction to the single-nucleon current. Although the latter cannot be derived by 
directly applying the Feynman rules, it can still be identified with the self-energy diagrams of Fig. 2 
(f)-(g). Thus in the next section we shall compute the renormalized one-body current using the 
new spinors and energy-momentum relation and then expand it in powers of the square of the pion- 
coupling constant /^/tti^. As we shall see, the unperturbed one-body current is thus recovered as 
the leading-order term and the first-order one coincides with the self-energy contribution. 
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3.3 Non-relativistic self-energy 

Before performing the expansion of the renormahzed current we briefly examine the non-relativistic 
limit of the self-energy diagrams in order to bring to light some differences that exist with respect 
to the fully relativistic case. 

The non-relativistic leading order of the self-energy current in eq. (10) is obtained by using the 
following prescriptions, which are valid in the static limit, 

E]^ -^^ m (47) 

75 ^ ~ o- • k (48) 

(49) 



1 1 



2 



J^2 _ jy^2 Jj.2 _|_ jy^ ^ 

Sf{P) ^ SrrriP) = ^ • (50) 

The electromagnetic form factor T^{Q) is also replaced by r^^(g), representing the usual non- 
relativistic one-body (OB) current acting over bi-spinors [g, Q. Using the above relations and 
performing the sums over spin and isospin indices, the s.e. current matrix element results 

j^,{p, h) C^ xl [^nr{p)Snr{P)r';,r{Q) + r{^.(Q)5nr(^)S„.(h)] Xh , (51) 

where Xp &iid Xh are two-components spinors. The non-relativistic self-energy function is given by 
The non-relativistic nucleon propagator in eq. (^) in the medium is then 

Snr{P0,P) = -2- + ^^nrip) ^ + --- = ^2 ^ , • (53) 

PO-J^ PO-^ PO-^ PO-^- 5]„r(p) 

As is well-known, this is a meromorphic function whose simple pole again defines the new energy 
of the nucleon in the medium, namely 

P^ 

enr(p) = :^ + 5^nr(p) , (54) 

since S„r(p) is a function only of p. 

Since the non-relativistic self-energy function in eq. (p^ ) does not depend on spin, the nucleon 
wave functions are not modified in the medium. In fact the corresponding Schrodinger equation in 
momentum space, including the self-energy, is simply given by 



|^ + E„„(p) 



4>nrip) =Po4>nrip) , (55) 



with the bi-spinor (pnrip) corresponding to the eigenvalue po = e„r(p)- 

The non-relativistic analysis of the nucleon self-energy current |l^ is much simpler than its rela- 
tivistic counterpart. Indeed, in the former the self-consistency is immediately achieved because the 
nucleon wave functions are not modified by the self-energy interaction and thus the first iteration of 
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Figure 5: Diagrammatical series for the one-body electromagnetic current with dressed external 
lines. 



the "Hartree-Fock" equations already provides the exact energy. Instead, in the relativistic frame- 
work the spin dependence of the self-energy modifies the Dirac-spinors, inducing an enhancement of 
the lower components. Moreover, the field-strength renormalization constant, namely the residue 
of the nucleon propagator in eq. (|53|) at the pole, in the non-relativistic case is just one. Hence the 
enhancement of the lower components and the spinors' field strength renormalization are genuine 
relativistic effects that are absent in a non-relativistic analysis where only the energy-momentum 
relation in the medium is altered by the self-energy diagrams. We shall show in next section that 
the two above-mentioned relativistic signatures can be recast as new pieces in the electromagnetic 
current acting over free spinors. 



4 Self-energy current to first order 

The particle-hole current matrix element in the HF approximation reads 



(56) 



where </>(p) are the new renormalized HF spinors introduced in eq. (|40|). Hence eq. ( pq ) represents 
the electromagnetic excitation of the p-h pair with dressed external lines corresponding to the sum 
of the diagrams shown in Fig. 5. In order to obtain a one-pion-exchange expression we expand 
eq. ( pq) in powers of the square of the pion coupling constant /^/m^ and single out the first-order 
term, that is, the one that is linear in f'^/m'^. We shall still refer to the current thus obtained that 
represents the OPE contribution as the "self-energy" current and, importantly, we shall show that 



it yields a finite contribution, free from the divergence problem of the current in eq. (11). 

To proceed we start by deriving the HF energy e(p) to first order in /^/m^. For this purpose 
we note that the functions A{P), B{P) and C{P) defined in eqs. (pl|-]23|) are of order 0{p /m^). 
Hence the following expansion of the Dirac mass in eq. ([29|) holds: 



fh{P) 



m 



l + AjP) 
1 - C{P) 



m [1 + A{P) + C{P)] + 0{f/mt) 



(57) 
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Inserting this into eq. (|35| ) for the energy and expanding again to first order in f /m^, we get 



Po = e(p) ^ [l-C{P) + B{P)]^p^ + m^ + 2m^[A{P) + C{P)] 
= Ep + AE{po,p) , 

where Ep = \/p'^ + m^ is the unperturbed free energy and 

A^(po, P) ^ ^{m^A + ElB - p^C) + Oif/mt) 
Ep 



(58) 



(59) 



is the first-order correction to the energy. Next we can insert the above value of po inside the 
argument of the functions A, B, C. Expanding the latter around the on-shell value po = Ep we get 



AiP) = Aipo, p) = AiEp + AE, p) = A{Ep, p) + 0{f/ml) ^ Ao{p) , 

where Aq(p) = A(Ep, p). Likewise to first order we obtain 

B{P) ^ B{Ep,p) ^ Boip) 
C{P) ^ C{Ep,p)^C^{p) . 



(60) 



(61) 
(62) 



Recalling that for P on-shell the functions A, B, C coincide, i.e., ^o(p) = Bq{p) = Co(p), we can 
insert these on-shell values into eq. ( p9[ ) and, neglecting terms of second order, i.e., 0{f'^/m'^), we 
finally arrive at the result 

PO = e(p) = Ep + ^Boip)im^ + i?^ - p^) + 0{f/mt) =Ep + ^2m^Bo{p) + 0{f/mt) . (63) 
^p ^p 

The above expression can be recast in terms of the on-shell value of the self-energy 

Eo(p) = 2mBo(p) , (64) 

which satisfies the relation 



S(^p,p)n(p) = So(p)n(p) 



(65) 



thus showing that the free spinors are eigenvectors of the self-energy matrix T,{Ep,p) corresponding 
to the eigenvalue So(p). Hence to first order in f'^/m'^, the HF energy in eq. ([6^) is found to read 



m 



e(p) ^Ep + ;^So(p) 



(66) 



in terms of the on-shell self-energy eigenvalue So(p). When compared with the non-relativistic HF 
energy in eq. (p^ ) it appears that, beyond the different expressions of the self-energy functions that 
hold in the relativistic and non-relativistic frameworks, an extra multiplicative factor m/Ep occurs 
in the relativistic case. 

Once the HF energy e(p) is known to first order in f'^/m^, we expand as well the renormalized 
spinors, namely 



lm[p) 



u(p,m(p)) 



Eip)+ m(p) 
\ 2Eip) 



X 
ap 



E{p)+m{p) 



X 



(67) 
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Actually, for later use in the calculation of the hadronic tensor, it is convenient to expand the 
spinor multiplied by the factor . / ™^P-' . 

V ^(p) 

Thus we start by expanding the Dirac mass in eq. ( [57| ) around the on-shell energy, obtaining 
m(p) = m [1 + Aoip) + Co(p)] + 0{f/ml) ~ m + So(p) , 



where use has been made of the on-shell self-energy in eq. (|64D. Likewise, using the HF equation 
(eq. (|35|)), the Dirac energy -E'(p) defined in eq. (|45|) is given by 



E{p) 



l~B 
l-C 



e(p)=.[l-i?o(p) + Co(p)] 



TTl 

Ep + — So(p) 
Ep 



m 



i?p + — So(p)~e(p). 



After some algebra the following first-order expressions are obtained 

£'n - m So 



'E + rfi 

2E 

1 



Im + Ep /^^-p 



2Ep 
1 



2Ep Ep 



i-S 



E + m m + Ep \ Ep 



Inserting eqs. ( |70| ) and (^) into the renormalized spinor in eq. (|67| ) we get 

X 



l m{p) 



n(p,m(p)) 



Im + Ep 
2E„ 



^ Ep-mT,o 
2£^p -Ep 



q-p 

m -|-£'ri 



i-I^U 



m So [rn m + Ep 



-Ep — m 



2E, 



X 



Ep + m cr • p 

2Ed m + E, 



-X 



Since 



(Ep -m)x 

- Ep + m /L X 
m + En 



(Ep7o - m) 



X 
o"P 

m + E, 



-X 



the first-order (in f /m^) renormalized spinor can be cast in the form 



^n(p,..(p)):.^/- 



/ N , So(p)Ep7o-m 
n(p) + ^^ ii— n(p) 



(69) 

(70) 
(71) 



(72) 



(73) 



We now also expand the field-strength renormalization function defined in eq. 



purpose we use eqs. (|68|j69|) , obtaining 

^2(P) 

which implies that 



^ ^ , , m^ dA ^ dB p^ dC 



-I Po=Ep 



Z2(p)~l + -a(p) 



(74) 
For this 

(75) 
(76) 
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with 



m 



^ dA 



dB 



2 dc 



l_ E 

Ep dpo ^ dpo Ep dpo 



a(p) = Bo{p) + 

Finally, collecting the above results and inserting them into eq. 

So Ep-fo 



(77) 



Po=Ep 



we get to first order 



l m{p) ^ 
Po{p)' 



IPJ 




u{p} + 



Er. 



2^ — ^i(p) + 2«(P)^(P) 



(78) 



Thus, within the OPE approach the renormalized HF spinors in the nuclear medium are charac- 
terized by two new elements with respect to the bare m(p). In appendix A we study in detail the 
second term on the right-hand side of eq. (|78t) . We show that {EpjQ — 'm)u{p) is directly connected 
with the negative-energy components in the wave function. Thus the second term in eq. ( [7^ ) in- 
volves -g^ ''2E "^ ^(p)i which can be written as a difference between two divergent terms, one being 
5i;'(p)S(p)u(p) and the other a kind of "renormalization" counterterm of the original self-energy 
current. The final result turns out, however, to be a finite quantity. 

Using the above expressions for the renormalized spinors, we now expand the renormalized 
one-body current matrix element to first order in /^/ttz^, getting 



/?TT.(p) /m(h) 



ii^F(P>h) 




u[p) 



PM _l_ PM 



So(h)Sh7o -m 



a(h) 



+ 



Eh 

So(p) Epjo - 



2Eh 



-p/i 



m 



y/EpEh 



Ep 2£', 
boB(P'h)+j^5i5(p,h)] . 



m^,, a(p)^„ 



n(h) 



(79) 



In eq. (79) the term Jq^ represents the usual one-body current matrix element evaluated with free 
spinors, i.e., 

j^B{p,h)=u{p)T^{Q)u{h) , (80) 

whereas jjig^ is a new renormalized self-energy current matrix element that includes the effects of 
the renormalization of the spinors. It can be decomposed according to 



JRSEiP^ h) = i^,5si(p, h) + JliSE2iP^ h) 



(81) 



where jjis^i embodies the correction arising from the new spinor solution of the Dirac equation in 
the medium and jjigE2 ^^^ °^^ stemming from the field-strength renormalization function \/^ in 
the medium. Their expressions are the following: 



JRSEliPM 
JRSE2iP^'^) 



u{p) 



Y^J. 



So(h) ^h7o 



a(h) -I- a(p) 



2Eh 

ioB(p» • 



_^ So(p) Ep-fo 



m 



Er: 



2Er, 



-p/i 



n(h) 



(82) 
(83) 



4.1 Gauge Invariance 

As shown above, the renormalized HF current matrix element, expanded to first order in f'^/m^, has 
been split into the usual one-body current and into a new renormalized self-energy current. In order 
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to be consistent with the one-pion-exchange model, we should add the contribution of the seagull, 
pion-in-flight and vertex correlation currents corresponding to the diagrams shown in Fig. 2(a-e). 
We point out once more that the self-energy diagrams (f) and (g), corresponding to insertions in 
external legs, should be disregarded in computing amplitudes (or currents) in perturbation theory. 
Instead, their contributions should be taken into account via renormalized energies and spinors as 
solutions of the relativistic HF equations. We have expressed the latter, to first order in /^/m^, in 
the form of a new current operator (denoted as RSE current). 
Then the total current in our model reads 

f{p, h) = j^^(p, h) + j^p^(p, h) , (84) 

where Jqpe embodies the seagull, pion-in-flight, vertex correlation and renormalized self-energy 
currents, namely 

3oPE = Js + Jp + Jv.c. + J RSE ■ v°5) 

In what follows we shall prove the gauge invariance of this current to first order in /^/m^. In 



so-doing it is crucial to take into account not only the full current in eqs. (84,85), but also the 
first-order correction to the energy of the particles and holes due to the self-energy interaction in 
eq. (pq). In other words, for a given momentum transfer q = p — h, the energy transfer should be 
computed as the difference between the particle and hole HF energies and not using the free values 
Ep and -Eh- Hence the energy transfer is 



m m 

— So(p) - -^- 



UJHF = Ep-Eh + — So(p) - — So(h) (86) 



and, in conformity, the associated four-momentum transfer is Q^p = {ujhf,^)- To make the 
following discussion clearer we denote with Qhf the HF four-momentum and with ujhf the HF 
energy transfer, to distinguish them from the on-shell values Q and uj. 

Divergence of the one-body current 

The divergence of the zeroth-order one-body current computed using the HF four-momentum 
transfer Qhf is given by 

Qhf,^j(^b(p, h) = u{p)QHF,^r'' (Qhf) u{h) = u{p)Fi (Qhf) 0HF«(h) , (87) 

where the nucleon vertex T^^[Qhf) is evaluated at the momentum transfer Qhf- Since Il(p)^u(h) = 
0, only the first-order contribution arising from the self-energy correction survives, namely 



QHF,MJos(p,h) = u(p)Fi(Q) 



^So(p) - T^So(h) 



7ou(h) 



In the above the Dirac form factor Fi is computed at the unperturbed value Q^, since we disregard 
second-order contributions. Note that the one-body current itself is not gauge invariant — its 
divergence yields a first-order term which turns out to be essential for the gauge invariance of the 
full current, as we shall see below. 

Divergence of the MEC 

As shown in |I|, after summing over the intermediate spin and isospin indices, the following 
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seagull and pion-in-flight Ip-lh matrix elements 



•u. , N /^ ^y. _, ^ V- "^ 2m{Q + 2H-2KY 



(89) 



(90) 

are obtained. These currents are already of first order in /^/rrt^; thus in computing their divergence 
one must use the unperturbed value of the energy transfer M, otherwise a second-order term arises. 
Then using the free Dirac equation and exploiting the kinematics we obtain (see Q for details) 



^ „. , N f ^v _. N V- "1 (2(K-P-m^) 2(K-H-m ^)] ,, , 

QMpM = -y^/l ^0,^Mv)ran E ^ I (p _ ^)2 _ ^2 - ( j^ _ ff )2 _ ^2 ] <^) 

(91) 
(92) 



Upon addition of these two equations the terms containing ^ cancel, leaving for the total divergence 
of the seagull and pion-in-flight the expression 

^/„ „. f ^v .. ,^ 2m i 2(K ■ H - m'^) 2(K ■ P - m"^) ] „, , , 



which can be further simplified by introducing the self-energy in eq. (15) for on-shell momenta, 
yielding finally 



Q/^iJ^ + J^) = 3^1 e3afen(p)r,rfe[S(p) - S(h)]«(h) . (94) 

Divergence of the vertex correlation current 

Starting from the Ip-lh matrix element of the v.c. current in eq. @ and applying the Dirac 
equation, we get (see Q for details) 

Q,jL.{P,h) = ^n(p)r,Fir,E^75(f- ^) (p ,^^2 _ ^2 75(f - ^Mh) 

(95) 

where we recognize again the expression of the self-energy matrix in eq. ( p!3| ) . Since the Dirac form 
factor can be split into an isoscalar and an isovector component according to 

F, = ^{Ff + F^Ts), (96) 
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which yields 



TaFlTa = 3Fi + iF^ esabTaU 



the divergence of the v.c. current written in terms of the self-energy function reads 



Qf^j^AP, h) = n(p) ( Fi + -F.^esabTaU ) P(h) - S(p)] u{h) . 



(97) 



(98) 



Comparing this result with eq. (94) we note that the e term above cancels with the MEC contri- 
bution. Hence 

Qf^ [JMEciP^ h) + i,^.,. (p, h)] = n(p)Fi [S(h) - S(p)] n(h) . (99) 

The above relation just expresses the Ward-Takahashi identity relating the full vertex correction, 
namely MEC plus v.c. (diagrams 2 (a)-(e)), to the self-energy matrix element. 



Divergence of the RSE current 

Here we obtain the divergence of the renormalized self-energy (RSE) current defined in eqs. (81,8^,83) 
For this purpose we first note that the divergence of J^se2 vanishes to first order because it is pro- 
portional to the OB current. Hence we write 



QtJ^JRSEiPM =u{p) 



T? /7i^o(h) E^-fQ - m So(p) -Ep7o -m 



E\i 2£'h 



2K 



Using the relation n(p) ^n(h) = and 



^i(p) ^7o^i(h) 



n(p)2(m7o - Fh)n(h) 
u(p)2(£'p - m7o)n(h) 



it is straightforward to obtain 



Q^i^5^(p, h) = n(p)Fi [S(p) - S(h)] tx(h) + n(p)Fi 






Sn(h) 



m 



So(p) 



u{h). (100) 



(101) 
(102) 



M(h) . (103) 



Remarkably the first term of this equation cancels with the divergence of the MEC plus the v.c. 
current, given by the Ward-Takahashi identity in eq. (p9[), whereas the second term cancels with 
the divergence of the OB current in eq. (p8|). We have thus proven that, within the present model 
up to first order in /^/m^, the total current in eq. ( |8^ ) satisfies the continuity equation, namely 

(104) 



QhfAJob + 3 MEC + JL. + Jrse) = • 



5 Relativistic self-energy responses 

In this section we derive the self-energy contribution to the nuclear response functions to first order 
in /^/m,^. We can thus compare with the results obtained in ||^] using the polarization propagator 
formalism. Our goal is to show that the results obtained in the two formalisms coincide, although 
they stem from different approaches, as we have previously emphasized. 
The one-body hadronic tensor in HE approximation reads 






^^ ^ ^-j^^(p,h)*j^p(p,h)5(^ + e(h)-e(p)) 



pSfi tptf^ 



{2TTr E{p)E{h)' 



(105) 
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where p = h + q and j^^(p,h) is the one-body HF current in eq. (pq) computed using the 
renormahzed HF spinors and HF energies of the particle and the hole. 

Next we use the expansions in eqs. ( |79| ) for the current jhf and (B6|) for the HF energies. In 
addition we expand the energy delta function to first order in /^/m^ according to 



5{u; + e(h) - e(p)) ~ 5(l^ + £^h - ^p) + 



d6{uj + Eh-Ep) 
duj 



^h ^p 



(106) 



Inserting all of these relations into the hadronic tensor and neglecting terms of second order we get 
for the diagonal elements of the hadronic tensor ^ 



W^^piu, q) ^ t^^^(w, q) + AW^'^sEi^, q) 



(107) 



(the summation convention is not in force in eq. (|107] )), where Wq'^{uj, q) is the usual OB hadronic 
tensor of a RFG, i.e., 



W^'B=T.T.j4^^i^\foBiPM'S{u^ + E^-E,) 



, , . (27r)3 EryEi, 
and AW^g^{u!,q) is the first-order self-energy correction 



(108) 



^WJ^'sE 



Sh 



m 



SS^ (2vr)3^pi?h 



&(P,h)|' 



23f? j^^ip, h)*j^5E(p, h)5{u; + E^- E^) 
d 



So(h) - — So(p) 



^h 



Er. 



dco 



5{uj + Ei,-Ep 



(109) 



In eq. ( |109| ) the first term corresponds to the interference between the OB and the RSE currents, 
while the second one, which shifts the allowed kinematical region because of the derivative of the 
energy delta function, is due to the modification of the nucleon energies in the medium. 
Carrying out the spin traces for the single-nucleon current 



EboB(p»l' 



1 



4m ^ 



Tr{r'^(Q)(^ + m)r'^(-Q)(f + m)} , 



(110) 



we get for the renormahzed self-energy response function 
d^h 1 



^WJ^'SE 



Tr 



(27r)3 AEpEh I 



|r^(Q) 



i:o(h) Ejijo - m a{h) 
Eh 2£'h 2 



+T^'iQ){P + m)r'^(-Q)(f + m) 



(^ + m)r^(-Q)(f + m) 
So(p) Ep-jo -m a(p) 



E, 



5^0 (h) ^h7o -m a(h) 



Eh 2£'h 

+ r'^(Q)(^ + m)r^(-Q) 

x5{uj + Eh-Ep) 



p 2£'p 

r''(-Q)(f + m 



So(p) Epjo -m a(p) 



Er, 



2K 



;f + m) 



^We only work out the diagonal elements of the hadronic tensor, since these are the ones that contribute to the 
unpolarized inclusive longitudinal and transverse response functions. 
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Figure 6: First-order contributions to the polarization propagator with self-energy insertions in the 
hole or particle lines. 



d^h 



1 



(27r)3 AE^Ey, 



Tx{T^'{Q){P + m)r^(-Q)(f + m)} 



m 



So(h) 



m 



So(p) 



d_ 
dto 



<5(cu + Ey,- Ep) . 



(Ill) 



More precisely, one should add two copies of eq. ( Ill ), one with the form factors appropriate to 
the proton and one to the neutron. 

In appendix B we show that this contribution to the response function is identical to the one 
obtained in [^ computing the imaginary part of the polarization propagator corresponding to 
the two diagrams of Fig. 6. This identity is not trivial: indeed in the case of the polarization 
propagator the response functions, with the Fock self-energy dressing the particle and the hole 
lines, are computed by representing the product of two nucleon propagators as the derivative of 
a single one to deal with the presence of a double pole in the integrand. In the present paper 
the problem has been solved differently. First the entire perturbative series with Fock self-energy 
insertions has been summed up and then the result has been expanded to first order, thus obtaining 
a finite first-order current operator. Because of the equivalence of these two procedures (as proven 
in appendix B) we are confident about the validity of the results we have obtained for the self-energy 
contribution to the nuclear responses. 

Finally we display the results obtained in our framework for the nuclear response functions^. In 
particular, we explore the impact on the latter of the new currents j^g^i ^^^ Jrse2 that arise from 
the enhancement of the lower components of the spinors and from the field strength renormalization 
^/Z2{pj, respectively. 

In Fig. 7 we show the on-shell self-energy (solid curve) and the field strength renormalization 
function (dashed curve) given by eqs. ( p^ ) and (^), respectively. The explicit expressions for 
T,q{p) and a{p) are derived in Appendix C. The T,q{p) obtained here is in good agreement with 
the results of ref. |13] and its effect on the single-particle energy in eq. ( |66| ) and on the effective 
mass in eq. (B3) is very small (less than ~ 3%). Note that a, which is linked to the current j^sE2 
of eq. (|83|), is much smaller than So(p)/£'p, which enters in jj^g^i through eq. (82). This implies 



^For clarity the results presented here have been calculated ignoring the hadronic [ttNN) form factors. 
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Figure 7: The on-shell self-energy T,Q{p)/Ep defined in eq. (64) (solid line) and the field-strength 
renormalization function a{p) given in eq. ([77D (dashed line) plotted versus the momentum p. 



that the effect of the enhancement of the lower components of the spinors dominates over the field- 
strength renormalization. This is very clearly seen in Figs. 8 and 9, where the various contributions 
to the longitudinal and transverse responses stemming from j^ig^i ^^'^ 3rse2 ^^^ displayed versus 
the transferred energy lo for momentum transfer q = 0.5, 1, 2 and 3 GeV/c. It is evident that the 
effect of j^5.^2 is negligible with respect to that oi j^^g^^. The separate contributions of the particle 
and hole self-energies are also shown: as q increases the contribution of the particle is suppressed, 
whereas the one for the hole survives. 

In Fig. 10 we compare the contribution to the longitudinal response of the renormalization of 
the wave functions (dashed) with that arising from the renormalization of the energies (solid). It 
clearly appears that the effect that is linked to modification of the energy due to the medium is 
the dominant one, the other being very small, especially for large values of q. Similar results are 
obtained for the transverse response. 



6 Conclusions 

As a more in-depth extension of the analysis carried out in [||] , in this paper we have again studied 
the set of one-pion-exchange operators that contribute to the electromagnetic responses of nuclei in 
the Ip-lh channel. Our main concern in this systematic approach has been to preserve fundamental 
symmetries such as Lorentz covariance and gauge invariance. As a zeroth-order approximation we 
have used the RFG, this model being motivated by its ability to handle issues of covariance and 
gauge invariance in a reasonably transparent manner. Its limitations are not expected to be too 



severe, since our entire focus is placed on the kinematic region of the QEP [14|. To go beyond 
the RFG implies a consideration of interactions between the nucleons, and as a first step, in this 
work we have limited the scope to one pion-exchange, postponing for future study a more complete 
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Figure 8: The contribution of the renormahzed self-energy current to the longitudinal response 
plotted versus uj. The nucleus is '^'^Ca with A;i?=237 MeV/c and the units are 10""*^ MeV~^ at g=0.5 
GeV/c (panel a), lO'^ Me\~^ at q=l GeV/c (panel b), lO'^ MeV"! at q=2 GeV/c (panel c), lO"'' 
MeV^"*^ at g=3 GeV/c (panel d). The separate contributions of the current jjisEi f°^ ^^^ particle 
(solid) and hole (short-dashed) and of the current Jjise2 ^^^ ^^^ particle (long-dashed) and hole 
(dotted) are displayed. 

treatment including heavier mesons. 

In assessing the role of the pious in the electromagnetic nuclear responses, the MEC are not the 
only contributions that arise in first-order perturbation theory. In fact the pionic correlations are 
intimately linked to MEC through the continuity equation and only when the full set of Feynmam 
diagrams with one pion-exchange is considered can one expect gauge invariance to be fulfilled. 

Concerning such correlations, in this paper we have paid special attention to the self-energy 
contribution. Indeed the iteration of the self-energy diagrams generates a "dressed" propagator 
in the medium. By the same token the self-energy generates "dressed" or "renormahzed" wave 
functions in the medium, solutions of an in-medium Dirac equation, where the self-energy plays 
the role of a mean relativistic potential. This equation also provides the dispersion relation linking 
the energy and momentum of the nucleon in the medium. Importantly, the new spinors should be 
multiplied by a renormalization function ^Z2(p). 

As the self-energy is generated by the interaction of a nucleon with the other nucleons in the 
medium, the solutions of the new Dirac equation should be used as input to re-compute the self- 
energy and so on. The exact answer is obtained through a self-consistent procedure. In this 
paper, however, we have just considered the first iteration: we have thus computed the self-energy 
current confining ourselves to first-order corrections to the energy and spinors - or, equivalently, 
to corrections linear in the self-energy - which correspond to diagrams with only one pionic line, 
in order to be consistent with MEC and the vertex correlation currents. 

Notably in the first-order expansion of the renormahzed spinors two new elements with respect 
to the non-relativistic approach emerge, one arising from the negative-energy components in the 
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Figure 9: Same as Fig. 8 for the transverse response. 

wave function produced by the interaction, and the second from the renormaUzation function \fZ2. 
These two corrections can be expressed as a new renormahzed self-energy current, j^g^^ acting 
over free spinors. Our results show that, for typical kinematics, the former constitutes a correction 
to the total self-energy of roughly 10-20%, whereas the latter is negligible. Moreover, while at low 
momentum transfers both particle and hole contributions in the former play a role, at high q only 
the hole contribution remains. 
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Appendix A. Properties of the spinor n(p,m(p)) 

The expansion 



^.(p,m(p))=.^- 



u{v) + 



T.qEjo 



m 



E 2E 



u{p) 



(112) 



transparently displays the effect of the self-energy on the free spinor u(p). Indeed the second term 
in the square brackets of eq. ( |112| ) corresponds to a negative-energy component with momentum 
p. In fact, the Dirac equation for a positive-energy spinor is given by 



(p • 7 -F m)M(p) = £;7ou(p). 



with E>d . 



(113) 
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Figure 10: The contributions of the first (sohd) and second (dashed) term in eq. ( |109D to the 
longitudinal response is plotted. The units and labelling are as in Fig. 8. 



Now if we apply the operator (p • 7 -|- vti) to the spinor (i?7o — r?T-)'u(p) we obtain 



(p • 7 -F m)(S7o - m)'u,(p) 



p • 7(i?7o - m)u(p) -I- m{Ejo - m)u(p) 
{-E'jo - m)p • 7u(p) -I- m{Ejo - m)u{p) 
(-£'70 - m){E^o - m)n(p) -|- m{Ejo - m)u(p) 
{—EjQ — m + m){E^o — rre)u(p) 
-£70(^70 -m)n(p). 



(114) 



Hence (£70 — m)u is an eigenvalue of the free Dirac Haniiltonian with eigenvalue —E. Therefore 
the operator Ejq — m transforms a positive-energy spinor u(p) into a negative-energy one. 

Moreover, it is useful to write down the spinor correction in an alternative form. Using the 
identity in eq. ( |114| ) we can write 

(115) 



(f - m)(£7o - m)ii(p) = 2^70 (£70 - m)n(p) . 
Dividing by 2E[f — m) we then obtain 



E'-fo — m 
2^ 



-HP) 



P — m 



70 (£70 -m)n(p) . 



(116) 



Hence the second term in the square brackets of the right-hand side of eq. ( |112| ) can be recast in 
the form 



So £70 - m 
E 2E 



n(p) 



1 



j^ TO --/o{Ejo - m)u{p) 
E p — m 

Sf{p) (1 - |7o) S(p)n(p) 



y-^ (1 - -70 ) S(p)n(p) 



(117) 
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The first term in eq. (|117| ), S F{p)T,{p)u{p) , corresponds to the one that enters in the original 
(divergent) self-energy current for a nucleon on-shell (eq. ([Tl|)). The subtracted term, with the 
factor ^7o inserted between the propagator and self-energy, cancels the divergence and yields a 
finite result. Thus it can be viewed as a "recipe" to renormalize the self-energy current. 

Appendix B. Renormalized self-energy response using the polar- 
ization propagator 

In ||l| we computed the first-order self-energy contribution to the polarization propagator corre- 
sponding to the two diagrams of Fig. 6. The latter splits into the sum of the two terms, with Fock 
self-energy insertions in the hole and particle lines respectively, and reads 



jjM^ 



-iTr 



Wf 



{T^'{Q)S{H)T.{H)S{H)T%-Q)S{H + Q) 



+r^{Q)S{H)V'{-Q)S{H + Q)J:{H + Q)S{H + Q)} , (118) 
where S{K) is the nucleon propagator in the medium, namely 



SiK) 



i^ + m) 
i^ + m) 



1 



J^2 _ j^2 _|_ ^g 

e{k - kp) 



+ 2Txie{kF - k)6iK^ - n?)e{ko] 



+ 



e{kF - k) 



J^2 _ j^2 _|_ ^g J^2 _ j^2 _ ^g^ 



(119) 



As in ||4[ the double poles appearing in the integrand of eq. (|118D are treated by writing the product 
of the two propagators as a derivative. Thus for the hole propagators one has 



S{H)i:{H)S{H) = 



P + m)Y.{H){JJ + m) 
d 



e{h - kp) 



+ 



e{kF - h) 



[H"^ — m? + ieY (H^ — fn"^ — iAo)"^ 



da 



a=0 



1 



H"^ — a — m? + ie 



+ 2me{kF - h)5{H'^ -a- m'^)e{ho] 



(120) 



and likewise for the propagators of the particle, S{P). This technique, however, cannot be applied 
to the divergent self-energy current in eq. (|l^), and accordingly we had to invoke renormalization 
techniques. 

Using the above prescription, subtracting the contribution coming from the vacuum and taking 
the imaginary part, the following expression was obtained for the self-energy contribution to the 
hadronic tensor Q: 



da 



+ 



dp 



d'h I^^{E^{a),h;Ep,p;q) 
d^h /oT(^h,h;S;(/3),p;g) 



^=0^ (2vr)3 4EhE'^{f3) 

where p = h + q, and modified energies for holes and particles, namely 



6iEi,{a) + qo- Ep)9{kF - h)e{p - kp) 
5{E^ + qo- E'p{/3))9{kF - h)e{p - kF) , 



K(«) 



yh? + a + m? 



p^ + (3 + rrfl 



(121) 

(122) 
(123) 
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have been introduced (in the above a and /? are real parameters). Finahy the functions /^^ are so 
defined 






Tr[r^(Q)(/f + m)^{H){P + m)T^{-Q){P + m)] 
Tr[r^(Q)(^ + m)T%-Q){P + m)S(P)(f + m)] . 



(124) 
(125) 



In order to prove the equivalence between the responses computed using the polarization propa- 
gator in eq. ( |121| ) and the result in eq. ( |111D , obtained using the renormalized current and energies, 
we perform the derivative with respect to a and fi. For a general function F{ho) we have 



dF{Ei^{a)) 



da 



a=0 



1 


\dF{ho)l 


2Ei, 


dho 



ho=Eh 



(126) 



Hence, interchanging the derivatives and the integral, we get for the hadronic tensor the expression 



W^" 



d^h 



1 



d 



+ 



+ 



(27r)3 4£^h-E^p dho 
d^h 1 



Ilo{ho,ti;Ep,p;q) 



2/in 



S{Eh + qo- Ep)e{kF - h)9{p - kp) 



ho=E^ 



(27r)3 AE^Ep 
d^h 1 d 



lfo"(i?h,h;i?p,p; 



_!_ d_ 



5{E^ + qo- Ep)9{kF - h)e{p - kp) 



(2^)3 AE^Ep dpo 
d^h 1 



Ioi{Eh,h;po,p; 



(27r)3 AEhE^ 



2po 
/o7(i?h,h;£;p,p; 



-I Po=Ep 

1 d 



2Ep dqo 



6{Ei, + go - Ep)eikF - h)e{p - kp) 
5{E^ + go - Ep)9{kF - h)e{p - kp) . 



In differentiating the function Iq^ in eq. ( |125| ), we first consider the term: 



d 



1 



dpo [2po 



(f + m)S(P)(/' + m) 
2m 



(127) 



(128) 



2K 



Er- 



Po=Ep 
(f + m) + 7o(f + m) + (Z' + m)7o 



Po=Ep 



+ 



^ ifI + m)^^iP + m) 



2Ey, 



dpo 



Po=Ep 



where use has been made of the results 



S(P)(f + m) 
;f + m)(f + m) 



So(p)(f + m) 
2m{P + m) , 



(129) 



(130) 
(131) 



which hold for P'^ on-shell and where So(p) is the eigenvalue of the self-energy for on-shell spinors. 
Next we should compute the derivative of the self-energy S(P). This function has the general 
structure given in eq. (|T^), and its derivative implies derivatives of the coefficients A, B, and C, 
namely 

d^(P) dA(P) dBlP) dCiP) 



dpo 



dpo 



dpo 



dpo 



-7 • p + B{P)jo 



(132) 



which must be evaluated for P^ on-shell. Using again eq. ( |13lD together with the identity 

{^ + m)j''{f' + m) = 2Pf'{P + m) (133) 
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we obtain, for P on-shell, 



2K 



1 ,„ ,(9S(P),^ 



1 
El 



dpo 
odAiP) ^^dB{P) odCiP) ^ „, ; 



dpo P dpi 







dpi 







Po=-Bp 



iP + m) 



a(p)(f + m), 



(134) 



where the definition of the function a(p) in eq. ( [77|) has been used. 

Finally, collecting the above results, the derivative in eq. ( |129| ) is found to read 



d 
dpo 



^(/' + m)S(P)(f + m) 



Po=Ep 



So(p) 

Er. 



joEp-m^^ , _\ , ^E, , _\7o-E^p-'^ 



2^. 



-(/' + m) + (f + m) 



2K 



+ a(p)(/' + m) 



(135) 



Hence the following expression 
'l^^{H,P,Q) 



d 
dpo 



2po 



TV r^(Q)(/r + m)r-(-Q) 



Po=Ep 



i;o(p)7o^p -m a{p 



Ep 2Ep 



+ 



[P + m) 



+ Tr:{rf'{Q){P + m)r''{-Q){P + m) 



So(p) -foEp -m a(p) 



Ep 2Ep 



(136) 



yields the derivative of Iqi{H, P, Q) with respect to po and the similar result 
'l^^{H,P,Q) 



d 
dho 



2ho 



J ho=Eh 



Tr r^(Q) 



i:o(h) -foE]^ - m a{h) 
Eh 2Eyi 2 



{P + m)T^{-Q){P + m] 



+ Tr r^(Q)(/f + m) 



So(h)7o£^h -m «(h) 
_ -Eh 2£'h 2 



r'^(-Q)(/' + m 

(137) 



holds for the derivative of I^q{H, P, Q) with respect to Hq. In addition, with the help of eq. ( |130D , 
we can write for on-shell momenta 

= 2mEo(h)TV{r'^(Q)(/f + m)r(-Q)(f + m)} (138) 



and, as well, 



/o^i"(/f,P,Q) = 2mEo(p)T¥{r^(Q)(/f + m)r(-Q)(/> + m)} 



(139) 



Finally the response functions are found as linear combinations of the diagonal components of the 
hadronic tensor, i.e. W^^. Using the above equations the latter reads 



lYf^p 



d^h 1 
(27r)3 AEhEp 



Tr r^(Q) 



£o(h) 7o£^h - m a(h) 
Eh 2Eh 2 



+T''{Q)iP + m) 



i:o(h) 7o£'h - m a(h) 
Eh 2Eh 2 



{P + m)T^{-Q){P + m) 
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+r'^(Q)(^ + m)r^(-Q) 



So(p) 70-Bp -m q(p) 



Ep 2Ep 



iP + m) 



+ r^(Q)(/f + m)r^(-Q)(/' + m) 



So(p) 70-Bp - m a(p) 



-E'p 2£'p 



+ 



Sh 



(27r)3 4£;h£;, 



x5(^h + qo- Ep)e{kF - h)e{p - kp) 
-Tr {r^(Q)(^ + m)T^'{-Q){P + m)} 



X So(h)^ - So(p)^ 1 ^5{E^ + go - Ep)e{kF - h)e{p - kp) , 
\ Eh EpI dqo 



(140) 



which coincides with the result in eq. ( |lllj ), obtained by computing the response functions using 
the renormahzed current and energy. 



Appendix C. On-shell self-energy and field strength renormaliza- 
tion function 



In this appendix we show in detail how to evaluate the on-shell self-energy in eq. (|6j) and the field 
strength renormalization function in eq. (^). They can be expressed in terms of the integrals I{P) 
and L^(P) in eqs. (|l3,|l8D as follows: 



2 , 



So(p) = 2mB{Ep, p) = -12m^ PoLo{po, p) - P-^abo, p) - m-^HPo, p) 

mi L 



Po=Ep 



(141) 



and 

a(p) 



Boip) + 



Er, 



2dA{po,p) „2dB{po,p) 29C{po,p) 

m 1- £/„ ;:; p 



-12m' 



mt 



dpo 

Lo{po,p) _ 
Po 



dpo 



dpo 



IiPo,p) + 



Po=Ep 

dLo{po,p) p dL3{po,p) m? dI{po,p) 



dpo 



Po dpo 



Po dpo 



Po=Ep 
(142) 



where we have used eqs. (Ell|2^) and the derivatives 



/ dA{po,p) \ 
\ dpo J 



Po=Ep 



f 



mt 



Po 



dLo{po,p) dL3{po,p) ^dI{po,p 



dpo 



dpo 



m 



dpo 



+Lo{po, p) - PoHpo, P)]pr,=E, 



PO=^p 



( dBipo,p) 

V dpo / p„=Ep 

/ dC{po,p) \ 

\ dpo Jpo=Ep 



f 



mt 



f 



mt 



Po- 



dLo{po,p) dL3{po,p) 29I{Po,p) 



dpo 



dpo 



m 



dpo 



Po=Ep 



Po 



dLo{po,p) dL3{pQ,p) ^dI{po,p 



dpo 
Po 



+Lo{po,p) -^3(P0,P) 

P 



dpo 



PQ=Ep 



m 



dpo 



(143) 
(144) 



(145) 
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By choosing the z-axis in the direction of p the angular integrals can be performed analytically, 
yielding 



/(^P,P) 

LsiEp^p) 
5I(po,p) 



1 f''P „ k , 7(n, k) + 2pk 
dk- — ;;- m ■ 



(27r)2 7o 4p^k l{p, k) - 2pk 

1 r'^F k 7(p, k) + 2pk 

(27r)2 Jo 4^ " 7(p, A;) - 2pA; 

1 /■'^^ J A;2 A;7(p, k) -f{p, k) + 2pA; " 

(2^ 70 1 2^ " 8p^E^ "" -f{p, k) - 2pk 



dpo 
dLo{po,p) 



1 r'^p k 

dk,^^{Ep-E^) 



dpo 
dL3{po,p) 



p^^E^ (27r)2 Jo 2pE^ 

1 f'^p k 

- 2 / dk—{Ep - El, 

^ h 2p 



1 



1 



dpo 



Po=Ep 
Po=Ep 



{2^f 



7(p, A;) + 2pk "f{p, k) — 2pk 
1 1 



1 [^F „ k 
dk 






7(p, k) + 2pk 7(p, k) — 2pk 
j{p, k) + 2pk 



(146) 
(147) 
(148) 
(149) 
(150) 



+7(p, A;) 
where we have defined the function 



(27r)2 Jo ip^Ek 

1 



7(p, k) — 2pk 



1 



7(p, k) + 2pk j{p, k) — 2pk 



7(p, k) = {Ep - E^f - np^ -k^-ml = 2rr? - m^ - 2EpEi^ . 
By replacing the above integrals in eqs. (141, 142|) we obtain 

1 + I!^ In 7(P^ k) + 2y>A; 



E„w = i!jC/"rf*l^ 



27r2m2 Jo Ey, 



Apk '~f{p, k) — 2pk 



and 



„(,).?!?:£ /"'dA-*' ^--^•' 



It is interesting to note that for large p- values the following limits hold 



lim a(p) = 



lim T,o{p) 

p—>oo 



^rnf f^^ „k^ ^rnffj^ , 2, Ep + kp 



(151) 
(152) 

(153) 
(154) 

(155) 
(156) 



where Ep = ykp + m^ ig the Fermi energy. For kp = 237 MeV/c the on-shell self-energy limit is 
~ 34 MeV. 
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